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•  A  generic  methodology  for  battery  cell  terminal  voltage  modeling  is  presented. 

•  Model  based  experiment  design  is  introduced  for  optimal  battery  cell  testing. 

•  The  optimal  excitation  signals  cover  the  entire  operating  ranges  of  non-road  HEVs. 

•  The  framework  is  purely  data-based  and  regards  relaxation  and  hysteresis  effects. 

•  Measurements  validate  the  methodology  applicable  to  different  cell  chemistries. 
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An  accurate  state  of  charge  (SoC)  estimation  of  a  traction  battery  in  hybrid  electric  non-road  vehicles, 
which  possess  higher  dynamics  and  power  densities  than  on-road  vehicles,  requires  a  precise  battery  cell 
terminal  voltage  model.  This  paper  presents  a  novel  methodology  for  non-linear  system  identification  of 
battery  cells  to  obtain  precise  battery  models.  The  methodology  comprises  the  architecture  of  local 
model  networks  (LMN)  and  optimal  model  based  design  of  experiments  (DoE).  Three  main  novelties  are 
proposed:  1)  Optimal  model  based  DoE,  which  aims  to  high  dynamically  excite  the  battery  cells  at  load 
ranges  frequently  used  in  operation.  2)  The  integration  of  corresponding  inputs  in  the  LMN  to  regard  the 
non-linearities  SoC,  relaxation,  hysteresis  as  well  as  temperature  effects.  3)  Enhancements  to  the  local 
linear  model  tree  (LOLIMOT)  construction  algorithm,  to  achieve  a  physical  appropriate  interpretation  of 
the  LMN.  The  framework  is  applicable  for  different  battery  cell  chemistries  and  different  temperatures, 
and  is  real  time  capable,  which  is  shown  on  an  industrial  PC.  The  accuracy  of  the  obtained  non-linear 
battery  model  is  demonstrated  on  cells  with  different  chemistries  and  temperatures.  The  results  show 
significant  improvement  due  to  optimal  experiment  design  and  integration  of  the  battery  non-linearities 
within  the  LMN  structure. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

A  novel  generic  methodology  for  non-linear  system  identifica¬ 
tion  and  optimal  model  based  design  of  experiments  (DoE)  of 
battery  cells  are  proposed  in  this  paper. 

The  control  strategy  of  hybrid  electric  vehicles  (HEV)  is  essen¬ 
tially  dependent  on  the  state  of  charge  (SoC)  of  the  used  traction 
battery.  The  state  of  charge  of  the  battery  is  not  measurable  on-line, 
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which  requires  an  estimate  of  the  actual  SoC  during  operation  [1], 
The  estimation  of  the  SoC  is  placed  in  the  battery  management 
system  (BMS)  and  is  often  only  based  on  the  open  circuit  voltage  of 
the  battery.  This  leads  to  big  estimation  errors,  since  the  non-linear 
behavior  of  the  battery  voltage  at  operation  is  not  regarded  with 
this  approach  [2],  Another  approach  is  to  integrate  the  battery 
current.  Disadvantageously,  current  offsets  are  accumulated,  which 
may  lead  to  estimation  errors  after  some  time.  The  third  possibility 
is  to  use  SoC  estimators  (e.g.  extended  Kalman  filter),  which  require 
a  model  of  the  battery  that  can  be  implemented  in  the  BMS  in  real 
time.  The  model  is  an  integral  part  of  the  BMS  and  describes  the 
non-linear  dynamic  behavior  of  the  battery  cell  terminal  voltage. 
The  SoC  estimation  accuracy  can  be  improved  only  if  a  precise  non¬ 
linear  battery  cell  model  is  used  in  the  SoC  estimator  [3], 
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Non-road  hybrid  electric  vehicles  and  machinery  (e.g.  con¬ 
struction  site  vehicles,  mining  vehicles,  ...),  compared  to  on-road 
vehicles,  usually  demand  higher  power  densities  and  load  dy¬ 
namics,  which  makes  the  modeling  of  the  batteries  more  compli¬ 
cated  [4],  The  non-linear  battery  effects  (e.g.  hysteresis,  relaxation, 
temperature  effects,  ...)  of  electrochemical  batteries  are  increased 
due  to  the  high  power  densities  [5], 

In  this  paper,  a  non-linear  data-based  battery  model  is  pro¬ 
posed,  which  can  be  used  for  the  purpose  of  accurate  SoC  estima¬ 
tion  in  non-road  vehicles.  Optimal  model  based  DoE  is  utilized  to 
optimize  the  excitation  signals  for  battery  measurements.  The 
optimal  excitation  signals  are  used  for  the  model  parameter  iden¬ 
tification  to  increase  the  accuracy  at  high  dynamic  demands. 
Furthermore,  due  to  the  experiment  design  and  data-based  struc¬ 
ture,  the  model  can  be  obtained  for  different  battery  cell  chemis¬ 
tries  within  a  reasonable  time  period. 

State-of-the-art  battery  models  and  DoE  are  reviewed  in  the 
following.  The  solution  approach  and  the  contributions  of  the  paper 
are  summarized  at  the  end  of  this  section. 

1.1.  State-of-the-art 

In  the  literature  three  main  types  of  battery  model  approaches 
are  mentioned: 

1.  Equivalent  circuit  models 

2.  Electrochemical  battery  models 

3.  Data-based  battery  models. 

Equivalent  circuit  models  (ECM),  as  depicted  in  Fig.  1  exem- 
plarily,  use  basic  electric  elements  in  order  to  model  a  battery  cell. 
The  main  intention  is  to  parameterize  the  model  using  physically 
interpretable  values.  Gao  et  al.  [5]  used  an  ECM  with  one  RC  circuit, 
that  accounts  for  non-linear  equilibrium  potentials,  rate-  and 
temperature-dependencies,  thermal  effects  and  response  to  tran¬ 
sient  power  demand.  Pattipati  et  al.  [6]  used  a  modified  equivalent 
circuit  model  for  SoC,  state-of-health  (SoH)  and  remaining  useful 
life  estimation  in  the  BMS.  The  high  power  density  application  in 
automotive  industry  requires  to  consider  the  behavior  of  the 
impedance  elements  (such  as  solution  resistance,  charge  transfer 
resistance,  and  Warburg  impedance)  in  a  simple  ECM  (see  Gomez 
et  al.  [7]). 

Electrochemical  battery  models  pursue  to  physically  model  the 
electrochemical  behavior  of  the  battery.  These  models  are  able  to 
simulate  the  chemical  states  of  a  battery  accurately  and  to  give 
insight  into  the  system  itself  [8],  while  the  disadvantage  is  that  they 
are  computationally  intensive.  Doyle  et  al.  [9]  modeled  a  lithium 
battery  cell  by  using  concentrated  solution  theory.  Partial 
differential-algebraic  equations  are  used  by  Klein  et  al.  [8]  for  state 
estimation.  Santhanagopalan  et  al.  [10]  used  a  single  particle  model 
(SPM)  for  SoC  estimation  with  an  extended  Kalman  filter.  However, 
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due  to  non-consideration  of  spatial  variation  of  the  states  in  the 
battery  cell,  the  success  of  the  SPM  model,  especially  at  high  cur¬ 
rents  or  long  duration  pulses,  might  not  be  valid  for  the  operating 
region  encountered  for  HEV  [11]. 

Data-based  models  are  a  useful  way  for  modeling  and  estima¬ 
tion  purposes,  although  in  general,  the  model  parameters  are  not 
physically  interpretable  [8],  Plett  [12,13,3],  used  a  data-based  non¬ 
linear  state  space  model  for  extended  Kalman  filter  SoC  estimation. 
The  model  takes  different  current  directions  into  account  and 
regards  a  “hysteresis  state”  as  well  as  the  relaxation  using  a  low 
pass  filter  on  the  current.  Battery  cell  chemistry  independence  of 
the  model  is  assumed.  Charkhgard  et  al.  [14]  applied  neural  net¬ 
works  to  battery  modeling.  Based  on  a  stochastic  fuzzy  neural 
network  [15,16,17],  Wang  et  al.  [18]  modeled  the  non-linear  dy¬ 
namics  of  current,  temperature  and  SoC  to  the  battery  voltage.  Xu 
et  al.  19]  used  the  same  model  for  SoC  estimation.  Hametner  et  al. 
[1]  applied  a  local  model  network  (LMN)  to  battery  modeling.  A 
LMN  is  composed  of  several  local  models  that  are  linear  in  their 
model  parameters  and  have  a  certain  area  of  validity  defined  by 
validity  functions  (see  e.g.  Refs.  [20,21,22]).  The  model  output  is 
non-linear,  due  to  the  non-linear  interpolation  of  the  local  linear 
models  (LLM).  The  LMN  is  constructed  by  an  iterative  algorithm, 
which  starts  with  one  global  linear  model  and  adds  a  LLM  to  the 
network  in  every  iteration  until  a  certain  threshold  is  reached 
(partitioning).  The  validity  of  the  new  LLM  lies  in  a  specific  form  in 
the  partition  space  of  the  model  and  depends  on  the  algorithm's 
strategy. 

All  of  the  mentioned  battery  models  require  measurements  to 
parameterize  the  model  parameters.  The  measurements  are  ob¬ 
tained  by  applying  a  current  excitation  signal  to  a  battery  cell  and 
recording  the  voltage  response.  Kroeze  et  al.  [23]  used  simple 
constant  discharge  and  charge  cycles  for  the  identification  of  an 
ECM,  while  Gao  et  al.  [5],  Chen  et  al.  [24]  and  Hentunen  et  al.  [25] 
made  use  of  a  discharge  pulse  excitation  signal  for  the  same  pur¬ 
pose.  A  discharge  pulse  excitation  signal  is  also  used  in  Ref.  [26  for 
a  model-based  estimation  of  an  electrochemical  battery  cell  model. 
More  advanced  ECM  (e.g.  linear  parameter  varying  models)  are 
identified  in  Refs.  [27,28,13]  using  a  pulse  profile,  regarding  charge 
and  discharge  mode.  Hu  et  al.  [29]  employs  an  asymmetrical  cur¬ 
rent  step  profile,  in  order  to  cover  a  wide  range  of  SoC  as  well  as  a 
wide  current  range.  This  profile  is  more  dynamic  compared  to  the 
other  excitation  signals.  An  example  for  low  dynamics  in  non-road 
applications  is  the  dynamic  Federal  Urban  Driving  Schedule  (FUDS), 
which  is  used  in  many  papers  as  validation  signal  (see  e.g.  Refs. 
[30,19,23,31]). 

Depending  on  the  model  approach,  the  design  of  the  experi¬ 
ment  plays  an  important  role  since  the  excitation  signal  has  a 
decisive  influence  on  the  parameter  estimation,  especially  for  data- 
based  model  approaches  [4,1  ].  Model  based  design  of  experiments 
can  be  used  to  create  optimal  excitation  signals:  The  goal  is  that, 
based  on  a  prior  model  of  the  process  (reference  model),  the  in¬ 
formation  obtained  from  measurements  is  maximized  and  pa¬ 
rameters  can  be  estimated  with  minimum  variance  [32],  In  this 
context,  the  Fisher  information  matrix  I  (FIM),  a  way  to  measure 
the  information  content  of  a  signal,  is  often  used  for  optimization  of 
an  excitation  signal.  Furthermore,  constraints  of  the  process  can  be 
regarded,  provided  that  the  reference  model  is  sufficiently  accurate. 

In  Ref.  [33],  a  local  model  network  based  generation  algorithm 
for  static  experiment  design  is  proposed.  Dynamic  experiment 
design  for  multilayer  perceptron  networks  by  choosing  optimal 
inputs  from  a  candidate  set  is  proposed  in  e.g.  Refs.  [34,35].  Stadl- 
bauer  et  al.  [36,37]  focused  on  the  dynamic  design  of  experiments 
based  on  multilayer  perceptron  networks.  Based  on  these  papers, 
Hametner  et  al.  [38]  proposed  a  design  approach  for  non-linear 
dynamic  experiments,  which  is  targeted  to  minimize  the  model 
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variance  of  dynamic  multilayer  perceptron  networks  as  well  as 
local  model  networks.  In  Ref.  [4]  the  benefit  of  model  based  design 
of  experiments  on  battery  modeling  is  investigated  and  compared 
with  a  dynamic  excitation  signal  from  the  literature.  Based  on  a 
linear  dynamic  model  and  predefined  current  levels,  in  Ref.  [1], 
model  based  design  of  experiments  is  proposed  for  optimal  SoC 
excitation  and  minimal  measurement  duration. 

1.2.  Solution  approach 

In  this  paper,  to  achieve  the  required  accuracy  of  the  battery  cell 
terminal  voltage  model  for  non-road  application,  the  data-based 
LMN  approach  is  used.  The  advantages  of  the  LMN  are  that  expert 
knowledge  can  be  included,  the  computational  effort  is  low,  no 
random  initialization  of  the  parameters  is  needed  and  the  LLMs  can 
be  interpreted  as  local  linearization  of  the  process  [39,40].  The  LMN 
is  constructed  using  the  local  linear  model  tree  (LOLIMOT)  algo¬ 
rithm,  while  the  LMN  structure  is  defined  by  corresponding  inputs, 
which  are  adapted  to  regard  the  battery  non-linearities  SoC, 
relaxation  and  hysteresis  as  well  as  temperature  and  current  in  the 
model.  Furthermore,  to  achieve  a  physically  appropriate  network, 
the  LOLIMOT  algorithm  is  adapted  to  start  with  a  prepartitioned 
network  and  is  prohibited  to  split  the  network  within  specified 
dimensions.  The  resulting  battery  cell  model  is,  due  to  less 
computational  complexity,  real  time  capable  and  can  be  obtained 
for  different  cell  chemistries. 

For  the  reduction  of  the  identified  model  parameter  variance, 
optimal  model  based  design  of  experiments  is  utilized  to  maximize 
the  information  content  of  measurements.  The  goal  is  to  sufficiently 
excite  the  system  dynamics  and  cover  the  entire  SoC  range  during 
measurement,  while  relaxation  and  hysteresis  effects  are  regarded 
additionally.  The  optimization  of  the  excitation  signal  is  based  on  a 
scalar  cost  function  J(l)  and  is  focused  on  real  load  ranges  that  are 
frequently  used  in  operation.  Since  battery  constraints  on  current, 
voltage  and  SoC  need  to  be  regarded,  a  constrained  optimization 
problem  results.  The  optimization  problem  is  solved  using  a 
gradient  based  algorithm. 

The  generic  methodology  of  non-linear  battery  model  identifi¬ 
cation  is  validated  by  means  of  real  measurements  on  different 
battery  cells.  The  cell  chemistry  independence  is  shown  on 
different  cell  chemistries  and  the  temperature  dependence  is  pre¬ 
sented  for  a  relevant  temperature  range  in  non-road  application 
between  12.5  °C  and  35  °C. 

In  summary  the  main  contributions  of  this  paper  are: 

•  A  generic  methodology  for  non-linear  system  identification 
applicable  for  different  battery  cell  chemistries  and  different 
temperatures. 

•  The  integration  of  corresponding  inputs  in  a  LMN  to  regard  non¬ 
linear  battery  effects  such  as  SoC,  relaxation,  hysteresis  as  well 
as  current  and  temperature  within  the  battery  cell  terminal 
voltage  model. 

•  Enhancements  of  the  LOLIMOT  construction  algorithm  to  use  a 
prepartitioned  network  instead  of  one  global  model  and  to 
prohibit  splitting  within  specified  partitioning  dimensions. 

•  Optimal  model  based  DoE  to  create  high  dynamic  excitation 
signals  to  significantly  increase  the  battery  model  accuracy  at 
high  dynamic  (non-road)  applications. 

•  A  real  load  cycle  analysis  for  identification  of  frequently  used 
load  ranges  in  operation  to  focus  the  excitation  signal  optimi¬ 
zation  at  these  ranges. 

The  remainder  of  this  paper  is  structured  as  follows:  In  Section  2 
the  LMN  approach  for  battery  modeling  is  described.  In  Section  3 
the  optimal  model  based  DoE  is  presented.  In  Section  4  the 


improvement  by  optimal  model  based  DoE,  the  comparison  of 
different  LMN  structures  and  the  dynamic  accuracy  of  the  proposed 
battery  model  are  demonstrated.  A  conclusion  is  made  in  Section  5. 

2.  Battery  modeling  using  local  model  network 

The  goal  of  this  section  is  to  develop  a  generic  methodology  for 
non-linear  identification  of  high  dynamic,  current— voltage  battery 
cell  models,  which  regard  non-linear  battery  effects  such  as  SoC, 
relaxation,  hysteresis  and  temperature  as  well  as  current  effects. 

The  remainder  of  this  section  is  as  follows:  First  the  general 
architecture  and  structure  of  LMN  is  described.  Then  the  con¬ 
struction  of  the  LMN  using  the  LOLIMOT  algorithm  is  explained  and 
finally  the  battery  model,  based  on  the  LMN  structure,  is  developed. 

2.3.  General  architecture  and  structure  of  LMN 

The  local  model  network  structure  is  basically  an  interpolation 
of  local  linear  models,  which  are  only  valid  in  a  certain  operation 
regime.  The  underlying  model  structure  of  the  LLM  is  chosen  to  be 
autoregressive  with  exogenous  input.  Based  on  the  global  non¬ 
linear  model  output  y,  the  global  parameter  vector  d  and  the 
input  variables  u ;,  the  regression  vector  q>  can  be  arranged  as 

<KM)  =  [yUr-.Mq  1]T, 

y=[y{k-\,0)...y(k-n,0)],  (1) 

u,  =  [ul(k-d)...ul(k  —  d  —  ml)],l=  1 

where  k  is  the  actual  time  instant,  n  is  the  output  order,  m;  is  the 
input  order  of  the  3-th  of  q  input  variables  and  d  is  the  dead  time. 
The  input  variables  u /  span  the  so-called  input  space  Q  of  the  model 
[41]  and  the  one  in  Eq.  (1)  regards  the  bias.  Note  that  the  local 
models  are  depicted  for  MISO  systems,  but  in  general  MIMO  sys¬ 
tems  are  possible  as  well  [40], 

The  optimal  model  order  of  the  inputs  and  outputs  can  be 
identified  by  different  methods:  The  theory  behind  goodness-of-fit 
(GOF)  techniques  as  well  as  recommendations  for  different  GOF 
techniques  are  given  in  Ref.  [42].  GOF  methods  from  the  statistical 
point  of  view  are  discussed  in  Ref.  [43],  A  trade-off  between  the 
goodness-of-fit  and  the  model  complexity  is  provided  by  Akaike’s 
information  criterion,  which  is  used  in  e.g.  Ref.  [44]  to  find  the 
optimal  order  of  an  ARX  structured  battery  model.  In  this  paper,  in 
order  to  obtain  the  optimal  model  order,  the  methodology  pro¬ 
posed  in  Ref.  [45]  is  used.  Different  selections  of  the  model  order 
are  analyzed  by  the  mean  squared  error  and  the  best  compromise 
between  complexity  and  accuracy  is  chosen.  The  obtained  model 
orders  are  given  in  the  results  section. 

The  output  of  the  i-th  of  M  LLMs  is  calculated  by 

yfk.  0)  =  </  (k,  t))dh  (2) 

where  6  =  [i?i...j>m]t  is  the  global  parameter  vector  of  all  M  local 
parameter  vectors  i?,.  The  weighted  aggregation  of  the  local  model 
outputs  forms  the  global  non-linear  LMN  output,  which  is  obtained 
by 

M 

y(M)  =  (3) 

i=l 

where  $,■  (/<)  is  the  validity  function  of  the  i-th  LLM.  The  validity 
function  is  calculated  by 
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Fig.  2.  LOLIMOT  scheme  for  two-dimensional  partition  space. 


m) 


w(M) 


(4) 


where  /x,-  (k,z)  is  a  Kernel  function  for  which  different  common 
functions  (e.g.  uniform,  triangle,  ...)  can  be  used  [40]  and  z  =  [Z\... 
z^]  includes  all  actual  partition  variables.  The  normalization  is 
required  in  order  to  obtain  the  sum  of  all  validity  functions  to  be 
equal  to  1.  The  f  partition  variables  in  z  span  the  so-called  partition 
space  Z  of  the  model  [22].  The  partition  space  Z  is  not  mandatory 
the  input  space  Q  [46],  though  input  and  partition  variables  may  be 
the  same. 


2.2.  Construction  of  LMN  using  LOLIMOT 


In  this  paper,  the  construction  of  the  local  model  network  is 
done  by  the  local  linear  model  tree  algorithm  presented  in  Ref.  [40], 
LOLIMOT  divides  the  partition  space  by  axis-orthogonal  splits  [41  ]. 
Iteratively,  the  worst  LLM  is  identified  by  the  quadratic  error  cri¬ 
terion  and  is  split  into  two  new  LLM  in  the  dimension  with  the  best 
improvement  for  the  global  model  output,  until  a  certain  threshold 
is  reached.  Here,  the  threshold  is  the  maximal  number  of  LLM  M. 

The  parameters  for  the  local  linear  models  are  obtained  by  using 
weighted  least  squares  (WLS).  In  order  to  avoid  numerical  effects, 
all  input,  partition  and  output  variables  are  normalized  from  0  to  1 
[40],  In  Fig.  2,  the  procedure  is  depicted  for  a  two-dimensional 
partition  space  with  more  detail:  In  every  iteration,  the  worst 
LLM  is  identified  and  is  then  split  in  all  possible  dimensions.  All 
divisions  are  checked  by  the  quadratic  evaluation  criterion  and  the 
best  alternative  is  chosen. 

In  the  LOLIMOT  algorithm,  the  Kernel  function  /q  (k,z)  is  chosen 
to  be  Gaussian  [40],  which  is  calculated  by 


w(M)  =  exp 


1  f(zi(k)-cn)2 

2  1  4 


fo(k)  ~  Ci*)2)) 

<  )r 


(5) 


In  case  of  axis-orthogonal  splits  of  the  partition  space  Z,  Cy  can 
be  interpreted  as  center  point  of  the  LLM  and  ay  is  the  individual 
standard  deviation.  The  standard  deviation  can  be  approximated  by 
ay  =  kffjAy,  where  Ay  is  the  spread  of  the  LLM  and  kaj  is  a  sharpness 
factor.  The  kernel  function  sharpness  kaj  is  defined  by  the  user  and 
can  be  interpreted  as  an  overlapping  factor  of  the  LLM  that 


influences  the  smoothness  of  the  non-linear  model  output.  Note 
that  the  optimal  sharpness  factor  depends  on  the  specific  applica¬ 
tion  and  is  different  for  the  partitioning  dimensions.  The  chosen 
sharpness  factors  k„j  are  given  in  the  results  section. 

Alternatively,  Jakubek  et  al.  [47  presented  an  iterative  con¬ 
struction  algorithm  for  LMN,  where  a  statistical  criteria  along  with 
regularization  is  used  in  order  to  allow  an  arbitrary  orientation  and 
extent  in  the  partition  space.  In  Ref.  [48],  a  residual  from  general¬ 
ized  total  least  squares  parameter  estimation  is  used  in  an 
expectation— maximization  algorithm  for  proper  partitioning  of  the 
LMN. 

Note  that,  LOLIMOT  has  less  implementation  complexity  due  to 
the  axis  orthogonal  orientation  compared  to  the  arbitrary  orien¬ 
tation  in  the  partition  space  of  the  mentioned  alternatives  and  is 
therefore  used  in  this  paper. 

2.3.  Battery  modeling  using  LMN 

In  this  section,  the  LMN  is  adapted  for  battery  modeling.  In  the 
following,  corresponding  inputs  to  regard  the  non-linear  battery 
effects  in  the  LMN  structure  and  the  LMN  structure  itself  are  pre¬ 
sented.  Furthermore,  the  enhancements  for  the  LOLIMOT  algorithm 
are  explained  to  increase  the  physical  meaning  of  the  battery 
model.  Fig.  3  depicts  the  flowchart  of  the  battery  model  construc¬ 
tion  process. 

2.3.1.  Corresponding  LMN  inputs  of  non-linear  battery  effects 

In  general,  distinction  must  be  made  between  physical  and 
chemical  non-linear  effects  in  electrochemical  batteries. 

In  Subplot  (a)  in  Fig.  4,  constant  charge  as  well  as  discharge 
curves  for  different  C-rates,1  measured  on  a  lithium  polymer  (Cell 
A)  and  lithium  iron  phosphate  (Cell  B)  chemistry,  are  depicted  over 
the  SoC.  The  SoC  is  obtained  by  integration  of  the  cell  current  [31  ] 
and  the  open  circuit  voltage  (OCV)  is  estimated  by  interpolation 
between  the  charge  and  discharge  curve  [49],  In  Subplot  (b)  a  cur¬ 
rent  step  sequence  at  different  temperatures  for  Cell  B  is  depicted. 

The  non-linear  physical  influence  of  current  and  temperature 
are  observable  in  Subplot  (b)  in  Fig.  4.  The  current  applied  to  the 
battery  cell  changes  the  battery  voltage  dynamically,  while  the 
temperature  affects  the  value  of  the  internal  resistance  of  the 


1  The  C-rate  is  the  quotient  of  current  [A]  and  battery  cell  capacity  [Ah], 
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Fig.  3.  Flowchart  of  battery  model  construction  process. 
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Fig.  4.  Subplot  (a):  Discharge  and  charge  curves  for  a  lithium  polymer  (Cell  A)  and  lithium  iron  phosphate  (Cell  B )  chemistry  at  different  currents.  The  open  circuit  voltage  (OCV)  is 
estimated  by  interpolation  between  the  charge  and  discharge  curve.  Subplot  (b):  Voltage  response  of  a  current  step  sequence,  similar  to  [29],  at  different  temperatures  for  Cell  B. 


battery  cell.  Due  to  lower/higher  temperature,  increased/decreased 
internal  resistance  leads  to  bigger/smaller  voltage  drops  if  current 
is  applied  to  the  battery  cell.  Hence,  the  current  as  well  as  tem¬ 
perature,  which  are  physically  measurable  on-line,  are  regarded  by 
the  corresponding  inputs  Ucurrent  and  zTernp,  respectively.  Note  that, 
temperature  gradients  are  small  compared  to  gradients  of  the 
current  and  therefore,  it  is  assumed  that  the  temperature  in¬ 
fluences  the  voltage  behavior  only  statically. 

From  Subplot  (a)  in  Fig.  4,  the  non-linear  influence  of  the  SoC 
on  the  battery  cell  voltage  is  observable.  Therefore  the  SoC  is 
regarded  by  the  corresponding  input  zs0c.  which  is  the  actual 
value  of  the  SoC.  The  SoC  can  be  obtained  by  integration  of  the 
current  or  alternatively  by  state  estimation  (see  e.g.  Refs.  [1,3]). 
Note  here,  that  the  SoC  can  only  be  changed  by  the  current  and 


thus  the  influence  on  the  system  behavior  is  assumed  to  be 
statical  only. 

The  dynamic  chemical  effect,  referred  to  relaxation,  is  not 
directly  measurable.  It  can  be  observed  in  Subplot  (b)  in  Fig.  4  at 
standby  current,  where  the  voltage  is  slightly  converging  to  a 
steady  state  value  [13],  The  time  axis  of  the  time  constants  for 
relaxation  and  current  effects  are  very  different.  Note  that,  relaxa¬ 
tion  also  acts  during  current  phases  [50],  Following  Ref.  13],  the 
relaxation  non-linearity  can  be  approximated  using  a  low  pass  filter 
on  the  current,  which  need  to  fulfill  certain  requirements: 

•  After  a  long  rest  period  and 

•  during  constant-current  discharge/charge,  the  relaxation  state 
needs  to  converge  to  zero. 
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Fig.  5.  Process  flowchart  of  the  DoE  to  obtain  the  optimal  excitation  signal  U.  Compare  blue  block  with  battery  model  identification  flowchart  in  Fig.  3.  (For  interpretation  of  the 
references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


The  relaxation  input  is  realized  by  the  dynamic  corresponding 
input 

^Relax  =  filt{^UQUTrent)  i  (6) 

where  Aucurrent  describes  the  change  rate  of  the  current,  in  order  to 
force  zero  DC-gain  of  the  filter  flit  ( • ).  In  this  paper,  the  filter  is 
designed  as  a  third  order  low  pass  filter  with  appropriate  chosen 
relaxation  time  constant  of  the  battery  cell.  The  time  constant  of  the 
relaxation  is  approximated  properly  on  the  basis  of  the  voltage 
converging  speed  at  standby  current  after  a  current  pulse  have 
been  applied  to  the  battery  cell. 

In  Fig.  4  in  Subplot  (a)  at  10%  SoC  and  in  Subplot  (b)  at  standby 
current,  different  shapes  of  the  charge  and  discharge  curve  are 
observable.  This  chemical  non-linear  effect  is  referred  to  as  hys¬ 
teresis  ]51].  A  first  order  hysteresis  input  is  used  as  corresponding 
input  ZHyst  for  the  hysteresis  effect.  The  value  is  obtained  by  the 
sign  ( • )  of  the  current  and  is  kept  at  the  last  value,  if  the  current  is 
zero. 

Other  known  chemical  effects,  e.g.  intercalation  effects,  con¬ 
ductivity  of  anode/electrolyte/cathode  or  concentration  gradients, 
which  are  not  directly  observable  in  Fig.  4,  are  neglected  due  to 
marginal  influence  compared  to  the  mentioned  effects. 


2.3.2.  LMN  structure  of  battery  models 

In  order  to  integrate  the  corresponding  inputs  to  the  LMN 
structure,  it  is  necessary  to  distinguish  between  static  and  dynamic 
influence  on  the  cell  voltage.  The  dynamic  corresponding  inputs 
need  to  be  included  in  the  input  space  Q  to  be  inputs  of  the  dy¬ 
namic  LLM.  The  static  corresponding  inputs  are  included  in  the 
partition  space  Z  to  be  used  for  partitioning.  Based  on  this  differ¬ 
entiation  and  the  mentioned  battery  effects,  the  partition  as  well  as 
input  space  can  be  denoted  by 

Z  =  [  Z\  Z2  Z3  ]  =  [  Zg0c  ZHyst  Zfemp  ]  -  r-j\ 

Q=[u  1  U2  U3]=[ZSoC  ^Current  ^Relaxl- 

In  order  to  regard  the  continuous  change  of  the  voltage 
depending  on  the  SoC,  the  static  SoC  input  zs0c  is  additionally 


regarded  in  the  input  space  Q.  Note  that,  for  observability  of  the  SoC 
using  a  Fuzzy-Observer,  the  SoC  is  also  required  to  be  an  input  of 
the  model  (c.f.  [1]). 

2.3.3.  User  defined  prepartitioning  of  the  LMN  structure  for  battery 
models 

The  choice  of  the  initial  partitions  of  Z  used  in  the  LOLIMOT 
algorithm,  influences  the  partitioning  decision  significantly.  In  or¬ 
der  to  prevent  physically  inappropriate  partitions,  the  LOLIMOT 
algorithm  is  forced  to  start  with  predefined  physically  appropriate 
partitions,  instead  of  one  global  partition.  Furthermore,  selected 
dimensions  of  the  partition  space  are  prohibited  to  be  split.  This  has 
the  advantage  of  reducing  the  computational  efforts  and  keeps  all 
partitions  physically  appropriate. 

The  hysteresis  input  ZHyst  refers  to  the  corresponding  cell  po¬ 
larization  and  is  therefore  split  into  charge  as  well  as  discharge 
mode.  Further  partitioning  within  this  dimension  is  prohibited. 
Based  on  Subplot  (b)  in  Fig.  4,  the  temperature  influence  within 
12.5  °C  and  35  °C  is  assumed  to  be  evenly  distributed.  Therefore,  the 
temperature  dimension  is  split  into  three  partitions  and  parti¬ 
tioning  is  also  prohibited,  because  a  simultaneous  split  of  charge 
and  discharge  mode  cannot  be  achieved  due  to  the  split  in  the 
hysteresis  dimension.  The  number  of  initial  partitions  is  therefore 
6,  while  the  partitioning  degree  of  freedom  is  limited  to  the  SoC 
dimension.  Alternatively,  expert  knowledge  can  be  used  to  set  up 
an  appropriate  prepartitioned  network. 

3.  Optimal  model  based  design  of  experiments 

The  goal  of  an  optimal  battery  excitation  signal  is,  that  the 
system  dynamics  are  sufficiently  excited  and  the  entire  SoC  range  is 
covered  during  measurement,  while  relaxation  and  hysteresis  ef¬ 
fects  are  regarded  additionally.  Batteries  are  tested  by  applying  a 
load  current  excitation  signal  U  to  a  battery  cell  and  the  cell  voltage 
is  recorded.  Usual  excitation  signals,  e.g.  Ref.  [29],  apply  interme¬ 
diate  current  steps  for  dynamic  excitation  of  the  battery  cell.  Non¬ 
road  applications  usually  require  higher  dynamic  excitation  signals 
to  sufficiently  include  high  dynamic  system  behavior  within  the 
measurements  [4], 
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For  that  purpose,  optimized  high  dynamic  excitation  sequences 
are  created  by  optimal  model  based  design  of  experiments.  Relax¬ 
ation,  hysteresis  as  well  as  constant  current  behavior  are  addi¬ 
tionally  regarded  by  constant  discharge/charge  current  with 
following  standby  current  in  front  of/after  the  optimized  high  dy¬ 
namic  excitation  sequence.  The  entire  SoC  operating  range  is 
covered  by  merging  optimized  excitation  sequences  for  evenly 
distributed  SoC  values  across  the  SoC  operating  range. 

In  the  following,  the  optimization  of  the  dynamic  excitation 
sequence  based  on  optimal  model  based  design  of  experiments  by 
simultaneously  regarding  the  battery  constraints  is  discussed.  Opti¬ 
mality  criteria,  based  on  the  Fisher  information  matrix  I,  are 
formulated  and  optimized  by  a  gradient  based  algorithm.  Load 
ranges  frequently  used  in  operation  are  identified  by  real  load  cycle 
analysis  and  are  regarded  in  the  optimization  in  order  to  improve  the 
obtained  model  quality  especially  for  these  ranges.  At  the  end,  the 
extensions  on  the  excitation  sequence  are  described  in  detail. 


3.1.  Formulation  and  calculation  of  optimization  criteria 


The  methodology  of  optimal  model  based  design  of  experiments 
is  used  to  create  excitation  signals  for  identification  of  models  with 
minimum  parameter  variance.  According  to  Ref.  [38],  an  excitation 
signal  U  can  be  optimized  using  the  FIM,  which  can  only  be 
improved  by  modification  of  the  system  inputs.  The  calculation  of  I 
is  based  on  the  partial  derivative  of  the  model  output  with  respect 
to  the  model  parameters.  For  the  LMN  battery  model  described  in 
the  last  section,  the  parameter  sensitivity  vector  ^  (/<)  can  be 
denoted  by 


m 


dy(k, ») 

dd 


(k)v(k,  e) 

<l>M(k)<p(k,  ft) 


k  =  1 , . . . ,  N , 


(8) 


where  M  denotes  the  number  of  LLM  and  (/<),  <p{k,0),  6  as  defined 
in  (4)  and  (1).  The  reference  model  can  be  obtained  by  two 
possibilities: 


•  A  LMN  model  is  available. 

•  A  different  model  (no  LMN  model)  is  available. 


In  this  paper,  a  linear  model  of  the  battery,  which  corresponds  to  a 
LMN  with  only  one  LLM,  is  identified  by  an  a  priori  available  excitation 
signal  and  is  further  used  in  the  following.  In  case  of  an  available 
complex  electrochemical  simulation  model,  simulation  results  can  be 
used  to  identify  a  LMN  as  a  reference  model,  alternatively.  Never¬ 
theless,  this  paper  should  show,  that  the  easy  approach  of  a  linear 
model  already  increase  the  model  quality  significantly. 

The  FIM,  which  is  based  on  \p{k),  can  be  defined  by 


I  =  dy(k.d)T 

n-2  Z—j 

k=  1 


00 


k=  l. 


■  N, 


(9) 


where  a  is  the  variance  of  the  measurement  noise  [52],  Using  a 
compact  notation  for  the  parameter  sensitivity  matrix 


W=[*T(1  )„./(N)]r,  (10) 

the  FIM  can  be  denoted  by 

1  =  \'VT'V.  (11) 

0-2 

For  the  optimization  of  the  excitation  sequence,  scalar  criteria, 
all  based  on  the  Fisher  information  matrix  1,  are  formulated.  Three 
common  criteria  exist  [52]: 
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A  -  optimality  :  JA  =  Jr{x  1  j  — >  min  (12) 

D  -  optimality  :  JD  =  det(X)  — >  max  (13) 

E  -  optimality  :  JE  =  Amin(X)  — >  max  (14) 

The  trace  of  the  inverse  of  the  Fisher  information  matrix  is 
subject  to  minimization  for  A-optimality,  D-optimality  is  targeted 
to  maximize  the  determinant  of  the  FIM  and  the  intention  of  the  E- 
optimality  is  to  maximize  the  smallest  Eigenvalue  of  the  FIM.  The 
advantage  of  D-optimality  is  the  higher  sensitivity  to  single 
parameter  covariances  compared  to  the  A-optimality  36], 

With  the  three  common  criteria  (12)— (14),  the  process  to  obtain 
the  optimal  excitation  signal  U  can  be  depicted  in  a  flowchart, 
which  is  shown  in  Fig.  5. 

3.2.  Formulation  of  the  optimization  problem 

The  main  goal  of  the  optimization  is  to  obtain  an  excitation 
sequence,  which  sufficiently  excite  the  battery  cell  by  high  dynamic 
currents,  while  the  battery  constraints  are  approached  simulta¬ 
neously.  In  general,  constraints  on  current,  battery  cell  voltage  and 
SoC  need  to  be  regarded: 

•  The  current  needs  to  be  within  minimum/maximum  values  in 
order  to  avoid  physical  damage  or  accelerated  life  time  reduc¬ 
tion  of  the  battery  cell. 

•  Battery  voltage  limits  of  battery  cell  chemistries,  such  as  lithium 
polymer  and  lithium  iron  phosphate,  must  be  met  in  order  to 
avoid  an  explosion  due  to  electrolyte  oxidation  or  fire. 

•  The  SoC  of  the  battery  cell  is  restricted  due  to  the  physically 
defined  capacity  of  the  battery  cell. 

In  this  paper,  two  possibilities  should  be  proposed  to  be  used  as 
current  constraints  in  the  optimization: 

1.  High  dynamic  excitation  between  physical  minimum  and 
maximum  current. 

2.  High  dynamic  excitation  between  load  (current)  ranges, 
frequently  used  in  operation. 

In  the  first  approach,  the  lower  and  upper  constraints  of  the  cur¬ 
rent  are  simply  defined  at  the  physical  minimum  and  maximum 
current.  In  the  second  approach,  which  is  used  in  this  paper,  the 
current  constraints  are  defined  by  analysis  of  real  load  cycles:  A  his¬ 
togram  with  a  defined  number  of  intervals  is  used  to  determine  the 
distribution  density  of  the  load.  The  interval  limits  of  the  histogram 
are  used  to  define  the  lower  and  the  upper  current  constraints,  while 
the  corresponding  distribution  densities  define  the  durations  within 
the  corresponding  constraint  ranges.  In  Fig.  6,  a  load  signal  depicted  in 
Subplot  (a)  with  the  corresponding  histogram  depicted  in  Subplot  (b), 
can  be  seen.  In  Subplot  (d)  the  current  constraints  are  depicted.  Note 
that,  the  advantage  of  the  second  approach  compared  to  the  first 
approach  is,  that  the  information  content  is  especially  increased  at 
operating  ranges  used  in  real  applications. 

The  voltage  constraints  are  included  indirectly  in  the  current 
constraints  by  limiting  the  minimal/maximal  currents  depending 
on  the  SoC,  since  the  reliability  of  the  output  constraints  to  accu¬ 
rately  regard  the  cell  voltage  limits  using  a  linear  reference  model  is 
low.  The  limitation  ratio  is  depicted  in  Fig.  6  in  Subplot  (c).  Note  that 
due  to  the  constraints  on  the  current  depending  on  the  SoC  and  the 
duration  within  the  current  limits,  the  maximal  deviation  of  the 
SoC  from  the  starting  SoC  value  is  limited.  Therefore,  an  adequately 
chosen  starting  SoC  value  defines  the  SoC  limits. 
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Fig.  6.  Construction  of  the  constraints  for  the  optimization.  Subplot  (a):  Scaled  real  load  cycle.  Subplot  (b):  Histogram  of  the  load  density  distribution  of  the  real  load  cycle.  Subplot 
(c):  Limitation  ratio  of  the  applied  minimal/maximal  physical  current  depending  on  the  SoC.  Subplot  (d):  Current  constraints  according  to  the  histogram  in  Subplot  (b).  Green:  upper 
current  constraints.  Red:  lower  current  constraints.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


Based  on  these  definitions  of  the  battery  constraints,  the  opti¬ 
mization  problem  can  be  formulated:  The  D-optimality  criterion  is 
very  much  used  within  optimization  problems,  due  to  the  invariance 
by  re-parametrization  of  the  model  as  well  as  the  higher  sensitivity 
to  single  parameter  covariances  compared  to  the  A-optimality  [32], 
Furthermore,  in  Ref.  [4],  significant  increase  in  batteiy  model  quality 
has  been  shown  by  using  the  D-optimality  criterion.  The  optimiza¬ 
tion  problem  is  therefore  stated  for  the  D-optimality: 

D-optimality:  max  det(X)  (15) 


The  second  term  in  Eq.  (16)  requires  to  calculate  the  single  de¬ 
rivatives  of  the  parameter  sensitivity  vectors  with  respect  to  the 
model  input,  which  furthermore  is  based  on  the  calculation  of  the 
derivative  of  the  regressor  q>(k,d)  as  defined  in  Eq.  (1)  with  respect 
to  the  input  U(r),  which  is  denoted  by 


d  <pT(k,0)  rdy(/<-l,0)  dy(k-n,0) 

dU(r)  [  dU(r)  dU(r) 

6ll* 6lk-\)r---5pl6(k-m,)r  0 


k>r, 


(18) 


<,  *.  f  Umin  (k)  <  U(k)  <  Um3X(k),  k=  1 , . . . ,  N 
S  t  \  (leRNxl 

where  U  =  ucurrent  corresponds  to  the  current  input,  which  is  the 
only  degree  of  freedom  of  the  optimization  problem,  since  the 
relaxation  input  and  the  SoC  are  directly  dependent  on  the  current. 


3.3.  Constrained  optimization 


In  this  paper,  the  gradient  descent  method,  which  has  low 
computational  efforts  and  sufficient  performance,  is  used  to  solve 
the  stated  optimization  problem.  The  gradient  g  =  [g(l)...g(N)]T  is 
composed  of  the  derivative  of  the  design  criterion  with  respect  to 
the  input  U(r )  for  all  observations  N.  The  single  gradients  are  ob¬ 
tained  by  the  trace  of  the  product  of  the  derivative  of  the  deter¬ 
minant  of  the  FIM  with  respect  to  the  parameter  sensitivity  matrix 
fit  and  the  derivative  of  fit  with  respect  to  U{r).  The  single  de¬ 
rivatives  are  computationally  intensive  [36  and  can  be  denoted  for 
the  r-th  observation  by 


d/oW  _  /d/D(W)  d»P  \ 

d U(r)  {  dWT  dU(r)J’ 


(16) 


where  8y  is  the  Kronecker  delta  function,  /  corresponds  to  the 
current  input  of  the  LMN  battery  model  and  the  recursive  calcu¬ 
lation  of  the  former  model  output  with  respect  to  the  input  U(r)  is 
given  by 


dy(M) 

dU(r) 


Sy{k,0)  dy(k-  l,fl) 

5y(k  —1,0)  d  U(r) 

s ^ 

recursive  calculation 

dy(k.d)  d  y(k-n.O)  dy(k,&) 

+  ay(fc-n,*)  dtj(r)  +Dt^rP 


k  >  r. 


(19) 


The  derivative  of  the  parameter  sensitivity  vector  with  respect 
to  the  model  input  can  then  be  stated  for  k  >  r  by 


dtfr(fc) 

dU(r) 


<PX  (Ic) 


d  <p(k,0) 
d  U(r) 


</>(k,  0) 


ddq  (k) 

hUirf 


<PM(k) 


d  <p(k,0) 
~dU[rT 


+ 


<p(k,  0) 


d<PM(k) 
d  U(r) 


Summarizing  Eq.  (20)  in  a  compact  notation 


(20) 


where  the  first  term  for  the  D-optimality  is  given  by  (cf.  [53]) 

^^=2JD(W)W[»PTW]_1.  (17) 

Due  to  the  appearance  of  the  inverse  of  the  FIM  in  Eq.  (17),  it  is 
important  to  be  mentioned,  that  the  FIM  is  required  to  be  a  regular 
matrix  with  full  rank. 


dq>  =  [d^r(l)  dxkT(N) 

dU(r)  [  d U(r)  ■■■  dU(r) 


(21) 


the  second  term  of  Eq.  (16)  is  obtained.  Using  the  calculated 
gradient  g  in  the  gradient  descent  method  with  an  adaptively 
adjusted  step  size  ?j,  the  update  of  the  excitation  sequence  can  be 
formulated  to 
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Enlargement  in  Figure  8 


Fig.  7.  Optimal  excitation  signal  obtained  from  optimal  model  based  DoE  for  Cell  B 
depicted  for  two  temperatures. 


u(”+1)  =  C/(")  +rrg{v),  v  =  0,1,2, ...  (22) 

s  t  !  umin  <  U("+1)  <  Umax 
\  l/e#xl 

The  excitation  sequence  is  updated  with  the  calculated  gradient, 
while  the  current  constraints  are  applied.  This  process  is  repeated 
until  no  further  improvement  is  obtained. 

Alternatively,  the  output  constraints  can  be  directly  included  in 
the  optimization  problem  formulation.  Several  approaches  are 
available  for  solving  the  resulting  non-linear  constraint 


Fig.  8.  Enlarged  view  of  the  optimal  excitation  signal  in  Fig.  7  for  Cell  B  at  two  tem¬ 
peratures.  In  the  second  subplot,  the  current  constraints  for  the  optimization  (green 
and  red  lines)  as  well  as  the  maximal  current  constraint  depending  on  the  SoC  (black 
dashed  lines)  are  included.  (For  interpretation  of  the  references  to  colour  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


optimization  problem:  Stadlbauer  et  al.  [36]  minimized  the 
(quadratic)  difference  between  the  gradient  of  the  design  criterion 
and  the  excitation  signal  increment  in  every  iteration,  while  the 
feasible  area  is  approached  simultaneously.  Hametner  et  al.  [38] 
realized  the  constrained  optimization  by  Lagrangian  multipliers. 
Other  approaches  to  be  applied  are  sequential  quadratic  pro¬ 
gramming  (see  Ref.  54])  or  numerical  multi-objective  optimization 
(see  Ref.  [55]). 

3.4.  Extensions  on  the  excitation  sequence 

As  introduced  at  the  beginning  of  this  section,  the  optimized 
excitation  sequence  is  extended:  In  front  of/after  the  sequence, 
constant  discharge/charge  current  with  following  standby  current 
is  added.  The  value  of  the  constant  current  is  the  maximum  current 
dependent  on  the  SoC  (see  Subplot  (c)  in  Fig.  6).  The  duration  of  the 
constant  current  is  related  to  the  longest  time  constant  of  the 
system.  The  duration  of  the  standby  current  need  to  be  a  multiple  of 
the  duration  of  the  constant  current  to  include  sufficient  informa¬ 
tion  about  relaxation:  Compared  to  Subplot  (b)  in  Fig.  4,  at  least  a 
duration  ratio  of  6  is  required  for  an  abated  voltage  degradation  at 
standby  current.  Therefore,  a  ratio  of  8  is  used  in  this  paper. 

The  final  excitation  signal  is  merged  by  extended  high  dynamic 
sequences.  The  sequences  are  created  for  evenly  distributed  SoC 
values  across  the  entire  SoC  range  of  interest.  Depending  on  the 
non-road  application,  the  SoC  range  of  interest  may  vary.  It  is 
important  to  mention,  that  SoC  deviations  from  the  desired  SoC, 
need  to  be  compensated  by  variation  of  the  durations  of  the  con¬ 
stant  currents. 

In  Fig.  7  the  obtained  optimal  excitation  signal  is  depicted  for 
Cell  B.  Subplot  two  shows  the  optimal  excitation  signal  with 
applied  current  constraints.  The  corresponding  voltage  response  is 
depicted  in  the  first  subplot  for  two  different  temperatures.  As  can 
be  seen  in  the  figure,  the  voltage  and  SoC  constraints  are  met.  Note 
that  the  voltage  as  well  as  SoC  constraints  are  indirectly  regarded 
due  to  the  current  constraints.  Also  the  entire  range  of  SoC  required 
for  non-road  application  is  covered. 

In  Fig.  8  an  enlarged  view  of  one  extended  excitation  sequence  is 
depicted.  In  the  second  subplot,  the  current  constraints  used  for  the 
optimization  as  well  as  the  maximal  current  depending  on  the  SoC 
are  included.  The  dynamic  current  within  the  lower  and  upper 
current  constraints  (green  and  red  lines)  has  been  optimized  by 
optimal  model  based  DoE.  In  case  of  a  non-linear  reference  model 
used  in  the  optimization,  the  spread  of  the  current  over  the  dy¬ 
namic  range  may  also  include  more  intermediate  steps  [4], 

4.  Results 

In  this  section,  the  results  of  the  generically  applicable  meth¬ 
odology,  proposed  in  this  paper,  are  discussed.  It  is  shown,  that  the 
framework  of  optimal  model  based  design  of  experiments  and  the 
LMN  approach  can  be  applied  to  different  battery  cells  with 
different  cell  chemistries.  The  improvements  reached  with  the 
optimal  excitation  signals  are  presented  (Section  4.2),  different 
LMN  structures  for  different  cell  chemistries  are  compared  to  show 
the  influence  of  the  non-linearities  (Section  4.3)  and  the  accuracy  of 
the  battery  model  T5  is  discussed  in  detail  on  different  figures 
(Section  4.4). 

For  reproducibility  and  comparability  of  the  measurements,  a 
defined  measurement  procedure  is  applied  to  the  three  examined 
battery  cells,  which  are  introduced  in  the  following:  Cell  A  is  a 
40  Ah,  prismatic  type,  energy  cell  using  the  lithium  polymer 
chemistry  with  a  nominal  voltage  of  3.7  V.  Cell  B  and  C  are  cylin¬ 
drical  type,  power  cells  using  the  lithium— iron-phosphate  cell 
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chemistry  with  a  nominal  voltage  of  3.3  V.  Cell  B  has  a  capacity  of 
4.4  Ah  and  Cell  C  1.1  Ah,  respectively. 

The  current  excitation  signals  are  applied  to  the  battery  cells 
using  a  battery  tester,  which  is  specified  in  a  voltage  range  of  -1  V 
to  10  V  and  a  current  of  ±240  A,  while  the  voltage  response  as  well 
as  the  cell  temperature  are  recorded.  A  Hoecherl&Hackl  source/ 
drain  module  is  used  for  closed-loop  current  control  and  a  National 
Instruments  USB  data  acquisition  board  is  used  for  the  measure¬ 
ments.  The  ambient  temperature  between  12.5  °C  and  35  °C 
(relevant  temperature  range  for  non-road  application),  which  lies 
within  the  temperature  limits  of  the  battery  cells,  is  controlled  by  a 
Votsch  climate  chamber.  The  measurements  are  executed  using  a 
LabView  control  software,  which  also  monitored  battery  states  and 
safety  issues.  The  complete  system  is  designed  for  up  to  10  kHz 
measurements  with  step  response  time  constants  of  the  source/ 
drain  module  for  less  than  200  (is  in  order  to  fulfill  the  high  dy¬ 
namic  requirements  of  non-road  applications.  All  measurements 
have  been  recorded  with  a  sampling  rate  of  100  Hz,  which  simul¬ 
taneously  is  the  sampling  rate  of  all  models. 

4.1.  Measurement  procedure 

In  order  to  guarantee  reproducibility  as  well  as  to  exclude 
different  initial  battery  conditions  as  the  reason  for  unexpected  and 
undefined  appearing  effects  in  measurements,  it  is  made  use  of  a 
defined  measurement  procedure.  All  different  cell  chemistries  are 
examined  based  on  this  procedure,  which  erases  all  short  term 
history  of  the  battery  cells  (cell  conditioning).  Furthermore,  the 
procedure  is  required  for  comparison  between  different  battery 
cells  and  chemistries,  respectively.  It  is  established  as  follows: 

1.  Initial  capacity  check  at  25  °C 

2.  Set  temperature  of  climate  chamber 

3.  Fully  charge  the  battery  cell 

4.  Discharge  until  initial  SoC  is  reached 

5.  Apply  excitation  signal 

6.  Repeat  3  to  5  until  all  excitation  signals  are  recorded 

7.  Repeat  2  to  6  until  all  temperatures  are  recorded 

The  initial  capacity  check  (see  enumeration  (1)),  is  applied  only 
once  to  obtain  the  true  capacity  and  in  order  to  reduce  the  influence 
of  the  short  term  history  on  the  battery  cell.  The  capacity  is 
calculated  by  the  last  of  8  applied  charge/discharge  cycles  with 
different  C-rates  from  4C  to  1C.  Additionally,  the  charge  efficiencies 
can  be  obtained  by  the  difference  between  the  quotient  of 
discharge  capacity  divided  by  charge  capacity. 

4.2.  Battery  model  quality  improvement  with  optimal  model  based 
design  of  experiments 

The  benefit  due  to  optimal  excitation  signals  should  be  shown 
by  comparison  of  two  local  model  networks  (models  I\  and  F2) 
with  the  same  configuration,  but  identified  using  training  data  from 
conventional  DoE  and  optimal  model  based  DoE,  respectively.  In 
order  to  exclusively  show  the  effect  of  the  training  data,  the  models 
are  only  parameterized  for  a  simplified  LMN  structure. 

The  partition  space  of  the  models  is  defined  by  Z\2  =  [zsoc  zHyst] 
since  the  temperature  is  hold  constant  at  22.5  °C  and  the  input 
space  is  defined  by  Q12  =  [^current]  because  SoC  and  the  relaxation 
parameter  are  neglected  in  the  models.  The  model  orders  of  current 
and  voltage  are  found  by  comparing  the  mean  squared  error  of  a 
linear  model  for  different  selections  of  the  order  [45].  A  good 
compromise  between  complexity  and  accuracy  is  chosen  to 

^Current  —  ^Voltage  —  :)- 


Table  1 

Summary  of  the  LMN  configuration  parameters  of  models  fi  and  f2,  which 
are  used  to  compare  the  benefits  of  optimal  model  based  DoE  to  conven¬ 
tional  DoE.  The  model  parameters  are  identified  by  different  training  data  at 
22.5  °C  temperature  (c.f.  Figs.  9  and  10). 


LMN  structure 

Configuration  parameters 

^Current  =  ^Voltage  =  5 

^12  —  [2SoC  2Hyst] 

k(j,SoC,12  =  0.75 

Q12  =  [^Current] 

^ff.Hyst  =  0.05 

M  =  30 

For  the  given  partition  space,  the  kernel  function  sharpness  of 
the  SoC  is  found  by  comparing  the  suitability  of  a  LMN  model  with 
10  local  linear  models  at  different  selections  of  the  sharpness  factor 
fcff.soc-  A  good  compromise  between  smoothness  and  strict  parti¬ 
tioning  is  chosen  at  k„: s0c,i2  =  0.75.  Note  here  that  a  smaller  value  of 
ka  reduces  the  overlapping  of  the  validity  functions.  The  sharpness 
factor  for  the  hysteresis  input  is  set  to  fc^Hyst  =  0.05  since  a  sharp 
separation  is  preferable,  because  smoothness  between  discharge 
and  charge  behavior  is  not  physically  reasonable.  The  number  of 
LLM  is  defined  to  the  threshold  M  =  30  since  this  is  the  real  time 
capability  limit  obtained  by  test  runs  on  a  battery  emulator  control 
unit2  (Beckhoff  Industrial  PC  C6515).  The  LMN  configuration  pa¬ 
rameters  for  models  F 1  and  F2  are  summarized  in  Table  1. 

Fig.  9  depicts  the  raw  data  of  the  training  data  obtained  by 
conventional  DoE  for  Cell  C  which  is  used  to  identify  the  parameters 
for  model  F The  step  profile  training  data  is  similar  to  the  asym¬ 
metrical  step  profile  mentioned  in  Ref.  [29]  and  is  more  dynamic 
compared  to  the  other  observed  excitation  signals  in  the  literature 
[4[.  Therefore  it  is  taken  for  comparison  with  the  optimal  model 
based  DoE  excitation  signal,  which  belongs  to  the  group  of  high 
dynamic  excitation  signals.  The  step  profile  is  constructed  by 
merging  alternating  intermediate  current  steps  in  charge  and 
discharge  direction.  Due  to  the  variation  of  the  step  durations,  the 
entire  SoC  range  is  covered.  Because  of  the  strict  DoE,  the  voltage 
behavior  is  not  directly  considered.  This  leads  to  voltage  violations, 
if  the  current  is  not  limited  by  expert  knowledge.  The  applied 
current  constraints  are  depicted  in  the  second  subplot.  Due  to  the 
limitation  of  the  current,  the  voltage  and  SoC  constraints  are  met. 
The  full  current  capability  of  the  battery  cell  cannot  be  achieved 
with  this  excitation  signal. 

The  parameters  of  model  F 2  are  identified  by  an  optimal  excita¬ 
tion  signal  which  is  obtained  from  optimal  model  based  DoE  for  Cell 
C.  In  Fig.  10,  the  raw  data  of  the  optimal  excitation  signal  is  depicted. 
In  general  only  current  constraints  are  applied  to  the  optimal  model 
based  DoE,  while  voltage  as  well  as  SoC  constraints  are  indirectly 
regarded  in  the  current  constraints.  The  training  result  of  model  r2  is 
included  in  the  first  subplot.  Compared  to  the  step  profile  in  Fig.  9, 
the  dynamics  of  the  optimal  excitation  signal  is  higher  and  the  cur¬ 
rent  in  the  middle  of  the  SoC  range  is  bigger,  though  all  constraints 
are  met.  Note  that  the  signal  is  obtained  by  adapting  the  method¬ 
ology  described  in  Section  3  for  the  capabilities  of  Cell  C. 

The  accuracies  of  the  models  are  validated  by  a  repeated  real 
(current-)  load  cycle,  which  represents  the  high  dynamic  re¬ 
quirements  of  non-road  applications.  During  the  load  cycle,  on 
average,  the  battery  cell  is  not  charged/discharged,  which  makes 
repetitions  of  the  cycle  possible.  The  entire  SoC  range  is  covered  by 
alternatingly  raising/lowering  the  mean  current  value  of  the  cycle, 
which  charges/discharges  the  battery  cell  on  average.  This  valida¬ 
tion  signal  is  referred  to  as  SoC  validation  signal.  Note  that,  the 


2  CPU  usage  of  the  battery  model  with  30  LLM  is  between  8  and  12%,  which  is  the 
critical  limit  for  real  time  usage  in  a  battery  emulator. 
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Fig.  9.  Step  profile  training  data  adapted  for  Cell  C  to  be  used  for  the  parameter 
identification  of  model  fi.  The  temperature  is  hold  at  22.5  °C.  The  training  result  of 
model  f  i  is  included  in  the  first  subplot. 


Fig.  10.  Optimal  excitation  signal  training  data  for  Cell  C  at  22.5  °C  to  be  used  for  the 
parameter  identification  of  model  f2.  The  training  result  of  the  model  is  included  in 
subplot  one. 

ambient  temperature  is  also  kept  constant  during  the  experiment 
(22.5  °C)  and  in  order  to  match  the  full  current  capability  of  the 
battery  cell,  the  real  load  cycle  is  scaled  to  the  maximum  current 
allowed  for  Cell  C.  The  simulation  result  of  the  SoC  validation  signal, 
is  depicted  in  Fig.  11.  The  SoC  trajectory  can  be  seen  in  the  last 
subplot. 

Comparing  the  envelope  curves  of  the  maximal  voltage  values  in 
Fig.  11,  it  can  be  noticed,  that  both  models  include  the  non-linear 
relationship  of  SoC  and  battery  cell  voltage,  though  the  average 
error  of  model  F i  (56.25  mV)  is  significantly  bigger  than  the  error  of 
model  F 2  (5.44  mV).  In  Fig.  12,  an  enlarged  view  of  one  real  load 
cycle  in  Fig.  11,  is  depicted. 

The  dynamic  behavior  of  the  models  I  j  and  F 2  are  similar,  as  can 
be  seen  in  the  first  subplot  of  Fig.  12,  but  model  F 1  has  a  large  offset 


Enlargement  in  Figure  12 
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Fig.  11.  Simulation  result  comparison  of  models  If  and  /2,  shown  on  a  SoC  validation 
signal  measured  on  Cell  C  at  22.5  °C.  The  signal  is  based  on  a  real  load  cycle,  which 
represents  the  high  dynamic  requirements  of  non-road  applications.  The  envelope 
curves  of  the  maximal  voltage  values  show  a  significant  higher  error  of  model  If. 

compared  to  model  TY  This  is  caused  by  the  lack  of  high  dynamic 
excitation  included  in  the  step  profile  training  data  of  model  F j.  The 
optimal  excitation  signal  leads  to  a  significant  increase  in  the  model 
accuracy  at  high  dynamic  excitation  of  the  battery  cells.  Note  that, 
though,  the  maximum  C-rate  of  the  signal  is  above  9C,  the  accuracy 
of  model  F2  is  still  satisfactory. 

In  the  following,  a  more  detailed  discussion  of  the  LMN 
approach  for  battery  modeling  is  based  on  battery  models  trained 
with  optimal  excitation  signals,  exclusively. 

4.3.  Comparison  of  battery  models  with  different  LMN  structures 
and  cell  chemistries 

In  this  section,  three  different  LMN  battery  models  (F3,  J4,  F 5) 
are  compared  with  each  other  based  on  two  different  validation 
signals.  The  models  differ  in  the  LMN  structure  to  show  the  influ¬ 
ence  of  hysteresis  and  relaxation  input  on  the  model  accuracy, 
while  the  temperature  is  regarded  simultaneously.  Furthermore, 
the  models  are  identified  for  two  different  battery  cells  (Cell  A  with 
lithium  polymer  and  Cell  B  with  lithium— iron-phosphate  chemis¬ 
try)  to  emphasize  the  cell  chemistry  independence  of  the  LMN 
approach. 

The  model  parameters  of  the  models  (F3,  F 4,  T5)  are  identified 
by  the  same  training  data  obtained  by  optimal  model  based  DoE 
(see  Fig.  7  for  an  optimal  excitation  signal  of  Cell  B  at  two  tem¬ 
peratures).  The  optimal  excitation  signal  is  measured  at  10  different 
temperatures  between  12.5  °C  and  35  °C  in  2.5  °C  steps  and  is  then 
merged  to  one  training  data  set  (raw  data  for  identification).  The 
same  is  done  for  Cell  A  and  B,  respectively. 

The  LMN  structures  and  configuration  parameters  of  the  models 
are  summarized  in  Table  2.  The  simplest  model  (T 3 )  does  not  regard 
the  hysteresis  and  relaxation  effects,  in  order  to  see  an  improve¬ 
ment  due  to  the  hysteresis  and  relaxation  input.  The  partition  and 
input  space  are  defined  by  Z3  =  [zSoC  Zcurrent  zTemp]  and 
Q3  =  [zSoC  “current]'  respectively.  Note  that  the  current  zcurrent  is 
only  interpolated  between  the  charge/discharge  behavior  at 
standby  current  and  does  not  refer  to  the  corresponding  polariza¬ 
tion  of  the  battery  cell.  The  initial  partition  of  the  current  input  is 
therefore  split  between  positive/negative  current  instead  of  charge/ 
discharge  mode.  Nevertheless,  in  model  F3  the  current  effects  are 
included  by  zcurrent-  The  SoC  zs0c  in  the  input  space  is  included  due 
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Fig.  12.  Enlarged  view  of  Fig.  11 :  The  comparison  of  models  f i  and  f2  at  one  real  load 
cycle  measured  on  Cell  C  at  22.5  °C.  Significantly  accuracy  increase  due  to  high  dy¬ 
namic  training  data  observable  for  model  f2. 


Table  2 

Summary  of  the  LMN  configuration  parameters  of  models  r3,  r4  and  T5.  The  influ¬ 
ence  of  the  hysteresis  as  well  as  relaxation  input  on  the  model  accuracy  is  illustrated 
by  different  Z  and  Q.  The  model  parameters  are  identified  by  the  same  training  data 
(c.f.  Fig.  7). 


LMN  structure  for  I3,  r 4  and  f 5 

r 3  =  [ZsoC  ^Current  zTemp]> 

^4  £4  =  [zSoC  zHyst  zTemp]> 

F 5  -Z5  =  [zSoC  zHyst  zTemp]» 

0-3  —  [zSoC  ^Current] 

Q-4  ~  [zSoC  ^Current] 

0-5  =  [zSoC  ^Current  uRelaxl 

Configuration  parameters  for  r3t  f4  and  f5 

msoc  =  1 

UReiax  =  3 

^Current  =  ^Voltage  =  5 

^-<t,SoC,345  —  ^ff.Temp  —  0*6 
kjr.Hyst  =  kff.Current  =  0.05 

M  =  30 

to  the  continuous  change  of  the  voltage  depending  on  the  SoC  as 
well  as  the  requirement  for  state  observer  design  (c.f.  [1  ]). 

In  Model  r4,  the  current  input  is  replaced  by  the  hysteresis  input 
in  order  to  regard  the  hysteresis  effect.  The  partition  space  is 
defined  by  £4  =  [zSoc  ZHyst  Ziemp].  while  the  input  space  is  Q4  =  Q3. 


Enlargement  (a)  in  Figure  14 


Note  that  the  complexity  of  the  model  does  not  change  in  this  case. 
Additionally  to  the  configuration  of  model  7’4,  the  relaxation  input 
is  included  in  the  input  space  of  model  f5.  The  input  space  follows 
to  Q5  =  [zSoC  ^current  UReiax].  while  the  partition  space  is  Z4  =  Z5. 
Note  here,  that  the  model  complexity  is  increased  by  the  relaxation 
parameter.  Model  F5  includes  all  relevant  electrochemical  effects 
discussed  in  the  manuscript. 

The  configuration  parameter  objectives  for  models  /’ 3,  T4  and  f5 
does  not  change  in  principle  to  those  of  models  T \  and  T 2  and 
therefore  mcurrent  =  ^voltage  =  5,  fcff.Hyst  =  0.05  and  M  =  30  can  be 
directly  taken  from  Table  1.  The  sharpness  factor  for  the  current 
input  is  set  to  current  =  fcff.Hyst  =  0.05,  since  a  sharp  separation 
between  charge  and  discharge  behavior  is  preferable.  The  optimal 
kernel  function  sharpness  of  the  SoC  and  temperature  is  found  to 
be  fc<7,SoC,345  =  kff.Temp  =  0.6,  because  the  temperature  input  slightly 
increase  the  smoothness  of  the  validity  function.  The  remaining 
model  order  of  the  SoC  input  is  set  to  ms0c  =  1.  due  to  no  dynamic 
influence  on  the  system  behavior  and  the  order  of  the  relaxation 
input  is  set  to  m Relax  =  3,  because  filt(  ■ )  is  designed  by  a  third  order 
low  pass  filter. 

The  validation  of  the  battery  models  is  done  with  a  temperature 
validation  signal  and  the  SoC  validation  signal :  In  the  temperature 
validation  signal,  the  aforementioned  real  load  cycle  is  repeated  at  a 
defined  SoC  level  until  the  ambient  temperature  is  heated  up  to  the 
upper  level  of  32  °C  and  cooled  back  to  the  lower  level  of  18  °C.  The 
same  procedure  is  applied  for  different  SoC  levels.  The  influence  of 
the  SoC  is  validated  with  a  slightly  changed  SoC  validation  signal  as 
introduced  in  Section  4.2:  In  order  to  strengthen  the  significance  of 
the  validation,  the  current  mean  value  of  the  cycle  is  more  often 
raised/lowered  in  between  the  signal.  Note  that,  the  average  charge 
direction  is  therefore  changed  more  often  compared  to  the  SoC 
validation  signal  used  in  Section  4.2.  Nevertheless,  the  entire  SoC 
range  is  covered.  The  temperature  validation  signal  and  the  SoC 
validation  signal  are  depicted  in  Figs.  13  and  15,  respectively.  A  more 
detailed  description  is  given  in  the  next  section. 

The  model  accuracy  is  evaluated  using  the  mean  squared  error 
(MSE)  and  a  normalized  root  mean  squared  error  (NRMSE)  as  given 
in  Eq.  (23) 


eNRMSE.  % 
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Fig.  13.  Resulting  plot  of  the  temperature  validation  signal  applied  to  Cell  B  and  the 
simulated  cell  voltage  using  model  r5.  The  real  load  cycle  is  repeated  at  different  SoC 
levels  while  the  ambient  temperature  is  heated  up  to  the  upper  level  of  32  °C  and 
cooled  back  to  the  lower  level  of  18  °C. 


where  Ns  is  the  number  of  samples,  y  is  the  measured  output  and  y 
is  the  simulated  output.  Note  that,  the  normalization  in  Eq.  (23)  is 
done  by  the  measured  output  values  of  the  specified  validation 
signal,  which  permits  the  comparison  of  the  different  battery 
model  accuracies.  In  Table  3  the  simulation  error  values  of  the 
different  LMN  architectures  and  cell  chemistries  are  depicted. 

In  Table  3,  it  is  observable  in  the  MSE  as  well  as  NRMSE  values, 
that  hysteresis  and  relaxation  input,  respectively,  increase  the 
battery  model  accuracy  of  both  cell  chemistries.  The  current  input 
in  model  T3  only  interpolates  between  the  charge/discharge 
behavior  at  standby  current,  which  is  not  physically  appropriate.  In 
contrast,  the  hysteresis  input  refers  to  the  corresponding  polari¬ 
zation  if  the  current  is  going  to  zero,  which  reflects  the  electro¬ 
chemical  behavior  of  battery  cells  more  than  the  current  input. 
Therefore  the  accuracy  of  model  f4  is  increased. 

Electrochemical  batteries  usually  have  very  different  time  con¬ 
stants  in  the  dynamic  behavior.  Measurements  need  to  be  recorded 
with  a  high  sampling  rate,  as  otherwise  the  fast  dynamic  behavior 
is  not  included  in  the  recorded  data.  Due  to  the  fast  sampling  rates, 
the  data  get  highly  correlated  and  therefore  numerical  problems 
due  to  ill-conditioning  arise  at  the  parameter  estimation  [56,57], 
This  phenomenon  is  referred  to  as  redundance  [56],  Since  the 
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Table  3 

Validation  error  values  for  the  LMN  battery  models  f3,  f4  and  f5  as  defined  in 
Table  2.  The  error  values  are  shown  for  Cell  A  and  B,  respectively.  Comparing  the 
models  with  each  other,  independent  from  the  cell  chemistry,  model  /  5  has  the  best 
model  accuracy,  model  f3  the  worst  and  model  /  4  denotes  intermediate  accuracy. 


Battery  cell 

Cell  A 

Cell  B 

Criterion 

MSE 

NRMSE 

MSE 

NRMSE 

Temp.  Val.  Sig. 

r3 

2.935e-4 

2.059% 

7.592e-4 

2.671% 

t4 

2.668e~4 

1.963% 

6.585e~4 

2.487% 

r  5 

2.447e~4 

1.752% 

5.236e~4 

2.218% 

SoC  Val.  Sig.  20  °C 

^3 

3.71  le-4 

2.537% 

5.652e-4 

2.933% 

t4 

3.559e~4 

2.485% 

4.8  lie-4 

2.706% 

T 5 

1.431e~4 

1.438% 

3.386e4 

2.270% 

SoC  Val.  Sig.  25  °C 

r  3 

- 

- 

4.622e~4 

2.928% 

C4 

- 

- 

4.094e-4 

2.756% 

t5 

- 

- 

3.109e~4 

2.402% 

SoC  Val.  Sig.  30  "C 

t3 

- 

- 

3.633e4 

2.847% 

r4 

- 

- 

3.393e4 

2.752% 

rs 

- 

- 

2.496e~4 

2.360% 

model  is  sampled  with  100  Hz,  it  is  very  difficult  to  identify  the 
corresponding  relaxation  time  constant  precisely  within  the  dy¬ 
namic  behavior.  Therefore,  further  improvement  is  achieved  by 
inserting  the  relaxation  input  to  the  input  space  of  model  F 5. 

Comparing  the  MSE/NRMSE  values  of  the  three  discussed 
models  {Fg,  J4,  F 5)  in  Table  3,  the  best  accuracy  is  obtained  by 
model  r5  while  the  worst  is  obtained  by  model  Fs. 

4.4.  Detailed  discussion  of  LMN  battery  model  r5 

In  the  following,  since  model  T5  includes  all  relevant  electro¬ 
chemical  effects  and  showed  the  best  accuracy  in  the  comparison, 
the  relaxation,  hysteresis  and  temperature  influence  should  be 
discussed  in  more  detail  based  on  model  f’5.  The  system  behaviors 
of  the  different  cell  chemistries  in  this  paper  are  similar  and 
therefore  the  resulting  plots  are  only  shown  for  the  measurements 
of  Cell  B.  In  Fig.  13,  the  temperature  validation  signal  is  depicted.  In 
subplot  three  and  four,  the  SoC  as  well  as  temperature  trajectory 
can  be  seen,  respectively. 


Enlargement  in  Figure  16 
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Fig.  15.  Resulting  plot  of  the  SoC  validation  signal  applied  to  Cell  B  at  three  different 
ambient  temperatures  (20  °C,  25  °C  and  30  ”C).  An  enlarged  view  of  the  marked  region 
is  depicted  in  Fig.  16. 


The  first  subplot  in  Fig.  13  shows  the  voltage  response  of  the  cell 
terminal  voltage  and  the  simulated  cell  voltage  using  model  Fs.  It 
can  be  observed,  that  the  model  output  changes  depending  on  the 
ambient  temperature.  The  model  mismatch  at  the  beginning  of  the 
dynamic  excitation  at  the  different  SoC  levels  comes  from  a 
changing  time  constant  of  the  relaxation  effect.  Since  the  time 
constant  of  the  filter  input  is  kept  constant,  this  effect  is  not 
regarded  in  the  model.  Nevertheless,  the  model  accuracy  is 
significantly  increased  by  the  relaxation  input  (c.f.  Table  3).  For  the 
SoC  transition  between  the  repeated  load  cycles,  low  constant 
current  is  used,  while  in  non-road  applications  usually  high  battery 
currents  (more  than  20C)  occur.  The  lack  of  sufficient  information  at 
low  constant  current  (below  2C)  within  the  training  data  explains 
the  small  model  mismatch  at  the  transition.  The  hysteresis  effect  is 
due  to  the  continuous  change  between  charge  and  discharge  not 


Enlargement  (a) 


2.8 


Enlargement  (b) 


Time  [min] 
Enlargement  (c) 


Time  [min] 


Fig.  14.  Resulting  enlargement  plots  of  the  temperature  validation  signal  (see  Fig.  13)  applied  to  Cell  B.  Enlargements  (b)  and  (c)  as  defined  in  (a). 
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Fig.  16.  Enlarged  view  of  marked  region  in  Fig.  15.  One  real  load  cycle  applied  to  Cell  B 
at  30  °C  ambient  temperature.  Simulation  result  of  battery  model  r5  at  applied  current 
rates  above  20C. 

directly  observable  in  Fig.  13.  Comparing  the  error  values  in  Table  3, 
the  improvement  due  to  the  hysteresis  input,  can  be  verified 
explicitly.  In  Fig.  14,  the  enlargement  of  the  marked  region  is 
depicted. 

In  enlargements  (b)  and  (c)  in  Fig.  14,  one  enlarged  load  cycle  is 
depicted  for  the  marked  regions  in  enlargement  (a).  It  can  be  seen, 
that  the  dynamics  required  for  non-road  applications  are  met  by 
the  model.  Comparing  enlargements  (b)  and  (c),  the  effect  of  the 
mentioned  changing  time  constant  of  the  relaxation  effect  can  be 
observed  by  the  slightly  bigger  voltage  error  in  enlargement  (b). 

The  validation  of  the  SoC  influence  is  depicted  in  Fig.  15.  The  SoC 
validation  signal  is  repeated  for  20  °C,  25  °C  and  30  °C  constant 
ambient  temperature  to  show  the  SoC  influence  at  different 
temperatures. 

The  invariant  error  on  the  changing  SoC  in  Fig.  15  Subplot  4, 
show  the  non-linear  consideration  of  the  SoC  by  the  battery  model. 
The  voltage  errors  at  the  different  temperatures,  depicted  in  Sub¬ 
plot  4,  are  all  within  a  small  tolerance  tube  of  less  than  a  maximum 
error  of  90  mV,  which  shows  the  good  accuracy  of  the  model  at 
different  temperatures.  Note  that,  the  current  rates  are  above  20C, 
which  can  be  seen  in  Fig.  16,  where  an  enlarged  view  of  one  load 
cycle  is  depicted  at  low  SoC. 

In  Fig.  16,  one  real  load  cycle  is  depicted  for  a  low  SoC  value  at 
30  °C  ambient  temperature.  Though,  the  voltage  behavior  at  low 
SoC  and  high  current  values  (above  20C)  is  strongly  non-linear,  the 
model  is  able  to  represent  the  underlying  behavior  of  the  battery 
cell,  accurately. 

5.  Conclusion 

In  this  paper,  a  generic  methodology  for  non-linear  system 
identification  of  electrochemical  battery  cells  is  presented.  The 
local  model  network  approach  is  applied  to  model  the  battery  cell 
terminal  voltage.  The  relaxation  and  hysteresis  effects  as  well  as  the 
SoC,  current  and  temperature  are  regarded  in  the  LMN  by  corre¬ 
sponding  inputs.  The  LMN  is  constructed  with  an  enhanced  LOLI- 
MOT  algorithm,  which  uses  predefined  partitions  and  reduced 
partitioning  dimensions  to  obtain  a  physically  appropriate 
network.  Optimal  model  based  design  of  experiments  is  utilized  to 
create  optimal  excitation  signals,  which  increase  the  information 
content  of  measurements.  The  excitation  signals  are  optimized  for 
load  ranges  frequently  used  in  real  load  cycles.  The  results  indicated 
that  a  simple  linear  reference  model,  used  in  the  optimization  of 
the  excitation  signal,  could  already  increase  the  model  accuracies 
significantly.  The  battery  model  errors  for  validation  signals  with 
currents  above  20C  have  been  below  3%  NRMSE.  The  high  model 
accuracies  could  be  achieved  for  different  battery  cell  chemistries 


and  temperatures.  A  defined  measurement  procedure  guaranteed 
the  reproducibility  as  well  as  comparability  of  the  results  and  the 
real  time  capability  of  the  battery  model  is  demonstrated  on  an 
industrial  PC. 
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